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ABSTRACT 

We present time-resolved optical spectroscopy of the dwarf nova CSS100603:112253- 
111037. Its optical spectrum is rich in helium, with broad, double-peaked emission lines 
produced in an accretion disc. We measure a line flux ratio Hel A5876/Ha= 1.49±0.04, 
a much higher ratio than is typically observed in dwarf novae. The orbital period, as 
derived from the radial velocity of the line wings, is 65.233 ± 0.015 minutes. In com- 
bination with the previously measured superhump period, this implies an extreme 
mass ratio of M2/M1 = 0.017 ± 0.004. The Ha and Hel A 6678 emission lines addi- 
tionally have a narrow central spike, as is often seen in the spectra of AM CVn type 
stars. Comparing its properties with CVs, AM CVn systems and hydrogen binaries 
below the CV period minimum, we argue that CSS100603:112253-111037 is the first 
compelling example of an AM CVn system forming via the evolved CV channel. 

With the addition of this system, evolved cataclysmic variables (CVs) now account 
for seven per cent of all known semi-detached white dwarf binaries with P orb < 76 min. 
Two recently discovered binaries may further increase this figure. Although the selec- 
tion bias of this sample is not yet well defined, these systems support the evolved CV 
model as a possible formation channel for ultracompact accreting binaries. The orbital 
periods of the three ultracompact hydrogen accreting binaries overlap with those of 
the long period AM CVn stars, but there are currently no known systems in the period 
range 67 — 76 minutes. 
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1 INTRODUCTION 

Cataclysmic Variable stars (CVs) typically consist of a white 
dwarf and a near-main sequence, late-type star in a close 
binary orbit, with the defining characteristic that mass is 
transferred from the companion to the white dwarf through 
Roche lobe overflow. (See |Warner| 1995 for a detailed review 
of cataclysmic variables.) 

A subset of CVs show quasi-periodic brightenings of 
several magnitudes, known as dwarf nova outbursts. These 
outbursts are thought to be the result of a thermal instabil- 
ity in the accretion disc which forms around the white dwarf 
(Osaki 1989), and detection through outbursts remains an 
important method through which new CVs are discovered 
(e.g. |Gans icke 2005; |Drake et al.||2009| ). Some dwarf novae 
also show superout bursts, which last longer than the normal 
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outbursts and are generally brighter. During superoutbursts, 
the accretion disc becomes axially asymmetric due to tidal 
interaction with the donor star. The resulting tidal stresses 
and stream interaction with the outer disc are observed as 
a photometric modulation in the light curve, known as su- 
perhumps. The superhump period is closely related to the 
mass ratio of the system and is a good proxy for the orbital 
period, being typically only a few per cent longer than the 
orbital period (e.g. |Patterson et al.||2005| |Kato et al.|2009 
|Gansicke et al.||2009| |Wood et al.|2011| ). 



Through the process of angular momentum loss, CVs 
evolve from long to short orbital periods, to a minimum pe- 
riod P^in- The period minimum occurs at the point where 
the donor star is driven so far out of thermal equilibrium by 
the mass loss that it can no longer shrink rapidly enough 
to attain thermal equilibrium, i.e. the thermal timescale of 
the donor star is longer than the mass loss timescale. Further 
mass loss therefore requires the orbit to expand to accommo- 
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date the donor, so the system evolves back to longer periods. 
(See e.g. [Rappaport et al.|1982||Kolb fe Bara»ffe|1999| ) If the 
donor is partially degenerate, further mass loss will cause its 
radius to grow, accelerating the process. The period mini- 
mum is observed to be P^in^ 80 minutes in hydrogen-rich 
systems ( Gansicke et al.||2009 ). 

By definition, the donor stars in CVs fill their Roche 
lobes, which for a given orbital period, defines the density 
of the donor star. For orbital periods much shorter than P^in 
the density of a Roche lobe filling star will be higher than 
for typical main sequence stars, which implies that the donor 
has to be evolved or degenerate. At present, only 42 ultra- 
compact accreting binaries with P or b< 76 min are known 
(Section [4|. Almost all of these belong to the small class 
known as AM CVn stars (see Solheim 2010 for a recent re- 
view). Their donor stars are generally deficient of hydrogen, 
allowing them to have very short orbital periods, observed 
to be in the range 5-65 minutes. The white dwarf primary in 
AM CVn systems accrete from either another, less massive, 
white dwarf ( |Paczy nski 1967; Faulk ner et al.||1972 ), from a 
low mass semi-degenerate helium star ( Savonije et al.| 1986 



Iben & Tutukov|[l987), or from an evolved main-sequence 



star which has lost its hydrogen envelope ( Podsiadlowski 
et al.|2003) . Only in HMCnc, the shortest period AM CVn 



system, has hydrogen been detected in the optical spectrum 
( Roelofs et al.|2010| ). There are three different AM CVn for- 
mation channels, characterised by the three types of donor 
star seen in these systems. Their relative importance has 
been a long-standing question. Population synthesis models 
suggest that the evolved CV channel is unimportant com- 
pared to the white dwarf and helium star channels, mainly 
because of the long evolutionary timescales required (e.g. 
Nelemans et al. 2004 2010 Yungelson et al., in prep.). These 



calculations are however sensitive to the specific form of the 
magnetic braking model used (|van der Sluys et al.|2005a|b| ) , 
and [Podsiadlowski et al. ( 2003| argue that a noticable frac- 



tion of AM CVn stars could form from CVs with evolved 
donor stars. 

In this paper, we present optical spectroscopy of 
CSS100603:112253-111037 (hereafter CSS1122-1110), a 
helium-dominated binary which is evolving along the CV 
track to become an AM CVn system. It was discovered in 
outburst on 2010 June 3 by the Catalina Real Time Tran- 
sient survey (CRTS, |Drake et al.p009) , at a magnitude of 
14.3. It was followed up photometric ally by VSNElQ ob- 
servers (Maehara 2010a b; Kato 2010), who reported a su- 
perhump period of P sh =65.434 ± 0.040 minutes (Kato et al. 
2010|). An optical spectrum of this system is available from 



the Sloan Digital Sky Survey (SDSS, |Abazajian et al.|2009| ), 
shown in Figure [I] It reveals a blue continuum dominated by 
Hei emission, along with obvious hydrogen Balmer emission 
lines. Subsequently C SS1122-1110 was a lso included in the 
SDSS CV catalogue of |Szkody et~aI1 ( |201l| p| who also noted 
the unusual helium-to-hydrogen line ratios in this binary. 



2 OBSERVATIONS 

2.1 Optical spectroscopy 

We carried out time-resolved spectroscopy of CSS 11 22- 11 10 
over three nights in 2011 March, using the FORS2 spectro- 
graph at the European Southern Observatory (ESO) Very 
Large Telescope (VLT) at Paranal, Chile. Observing condi- 
tions were very good, with sub-arcsecond seeing throughout. 
We used the GRIS1200R+93 grism with a I" slit and 2x2 
binning of the 2kx4k MIT CCD. The resulting spectra cover 
the spectral range 5872 - 7370A with an average dispersion 
of 0.73A per binned pixel and a full width at half maximum 
(FWHM) resolution of 1.7A. 

We obtained 52 spectra in total over the three nights, 
covering approximately two binary orbits on the first night, 
two on the second night and one orbit on the last night. 
All spectra were taken using an exposure time of 420 sec- 
onds, except for one exposure on the second night which was 
aborted after 177 seconds due to a temporary guiding error 
of the telescope. The spectra were reduced using the STAR- 
lini^] packages KAPPA and FI GARO, and optimally extracted 
using PAMEL^] |Marsh[^989| ) . 

The wavelength calibration is based on helium-argon- 
neon arc lamp exposures which were taken during the day. 
We identified 25 reliable arclines over the spectral range and 
fitted them with a a fifth order polynomial. The root mean 
square (rms) of the fit residuals was 0.012A. In order to 
account for a small amount of flexure of the instrument, we 
shifted the spectra so that the strong night-sky emission line 
at 6300. 304A appears at its known wavelength. The typical 
shifts required were 0.1 — 0.2 pixels. 

The instrumental response and flux calibration were 
calculated from a spectroscopic standard star observed at 
the start of each night as part of the ESO calibration pro- 
gramme. The spectra from the first and third nights were 
calibrated using HD49798 and the second night's spectra 
were calibrated against Hiltner600. The spectra are not cor- 
rected for telluric absorption or slit losses, so the flux cali- 
bration is not absolute. However, the continuum slope and 
the relative calibration over the three nights are reliable. 



2.2 Photometry 

The Catalina Sky Survey (CSS) has visited the field contain- 
ing CSS1122-1110 256 times since the start of the survey. 
Our target was undetected on 82 of these images, and we 
rejected a further six for which the measured photometric 
error exceeded eight per cent (compared to the mean error 
of 2.3 per cent). This left 168 reliable photometric measure- 
ments. The median time between CRTS visits is 19 days, 
with typically 2-4 observations per visit, taken over 30- 
40 minutes. The light curve, spanning 8.1 years, is shown 
in Figure [2] Only one outburst was seen during this time. 
We measure an average unflltered quiescent magnitude of 
20.2 ± 0.5, in good agreement with the SDSS magnitudes. 
We also carried out differential photometry on our 



3 http:/ /starlink.jach. hawaii.edu/starlink 

1 http://www.kusastro.kyoto-u.ac.jp/vsnet/ 4 PAMELA is included in the STARLINK distribution 'Hawaiki' and 

2 Referred to as SDSS J112253. 3-111037.6 later releases. 
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Figure 1. Main panel: Optical spectrum of CSS100603:112253-111037 from the SDSS, with the strongest emission lines labelled. The 
SDSS spectrum was smoothed by a 15 point boxcar for display purposes. The gray line overplotted on the spectrum is a blackbody of 
T = 10 400 K, fitted to the SDSS and GALEX fluxes (red dots). The offset between the photometric and spectroscopic measurements 
is likely due to the variability of the source. The spectrum and the photometric measurements were corrected for Galactic extinction. 
Inset: Average spectrum from our VLT/FORS2 observations. 



15 - 
-S 16 - 

§ 18- 

J 19 " 
§20-1 

21 - 

22 - 



I I f f 1 



Table 1. Spectral energy distribution. 
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Figure 2. CSS lightcurve spanning eight years of observations. 
The average quiescent magnitude is 20.2 db 0.5 and the outburst 
reached a magnitude of 14.3. 



VLT/FORS2 acquisition images, using iraiQ The acquisi- 
tion was done using the RJ3PECIAL+76 filter, a broadband 
filter centred on 6550A with a FWHM of 1650A. It there- 
fore includes all four lines in our spectrum. The photometry 
showed that the object brightened by 0.15±0.02 magnitudes 
over the three nights of observations. 
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3 ANALYSIS AND RESULTS 
3.1 Spectral energy distribution 

The ultraviolet (UV) and optical magnitudes of CSS 1122- 
1110, as given by the GALEX (Galaxy Evolution Explorer, 
Mart in et al.|2005 ) Data Release 6 and SDSS Data Release 



8 ( |Aihara et al.||2011 



are shown in Table [T] The system 
is undetected in the GALEX far-UV band. Correcting for 
Galactic extinction towards the source, E(B — V) = 0.048 



(Schlegel et al. 1998), the photometric measurements are 
well described by a T = 10400 ± 300 K blackbody spec- 
trum, shown as a gray line in Figure [I] The flux difference 
between the photometric and spectroscopic SDSS measure- 
ments is due to the variability of the source. A lower limit on 
the distance to the system can be set by assuming that the 
continuum emission originates from the white dwarf only, 
with no disc contribution. Assuming a typical white dwarf 
radius of O.Oli?©, this gives d ^ 345 pc. 



5 IRAF is distributed by the National Optical Astronomy Obser- 
vatory (NOAO), which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under cooperative 
agreement with the National Science Foundation. 



3.2 Emission line profiles 

We show the average line profiles of the four emission lines 
visible in our FORS2 spectra in Figure [3] Unfortunately, 
much of the Hei A 5876 emission line is lost from the blue 
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end of the CCD, so this line is not usable for analysis. The 
lines are broad and display a classic double-peaked profile, 
revealing the presence of an accretion disc. The Ha and 
Hei A6678 lines additionally have a narrow emission compo- 
nent at the line centre. The Hei A7065 and Hei A7281 lines 
show no such component. Such a 'central spike' is never 
seen in the spectra of CVs, but many A M C Vn systems dis - 
play this feature^ e.g. GP Com ||Smak||l975| |Marsh||l999| ), 
V396Hya (CE315, |Ruiz et al.||2001| ) "and SDSSJ124058.03- 
015919.2 ( |Roelofs et al.||2005| ). In these systems the central 
spikes display radial velocity variations consistent with emis- 
sion from near the white dwarf ( |Morales-Rueda et al.|2 0Q3 ) . 
A similar narrow emission line originating from near the 
white dwarf is seen in some detached post-common-envelope 
binaries (e.g. LTT560, Tappert et al. 2011b but see also 



Table 2. Line equivalent width variation over the three nights 
of observations. 



Tappert et al.||2011a| and Fig. 8 of IGansicke et al.||2004). 



The white dwarf in LTT 560 is accreting from the wind of 
its companion, and the emission is thought to originate in a 
chromosphere or corona around the white dwarf which forms 
as a result of the accretion activity. 

The central spike in CSS 1122- 11 10 is too weak and its 
variations too small to reliably measure the radial velocity 
amplitude with a moving Gaussian model. The trailed spec- 
tra in Figure [3] show that it is stationary at the line centre 
to within one pixel (33 kms -1 ). Future observations will 
greatly benefit from using a shorter exposure time per spec- 
trum to reduce the amount of orbital smearing of the spike, 
which will allow us to track the small amplitude variations 
with greater precision. 

We observed an increase in line and continuum flux over 
the three nights of observations. In order to quantify the 
variability, we shifted the emission lines to their zero velocity 
position and calculated the average of each night's spectra. 
The equivalent widths, as measured from these average spec- 
tra, are shown in Table [2] The equivalent widths increase by 
a factor > 1.5 from the first to the third night. The line flux 
increased by a factor of 1.6 during the same time, and the 
continuum, as measured from a line- free part of the spec- 
trum (50A centered on 6 100 A), brightened by a factor of 
1.2. The brightness increase is consistent with the magnitude 



change measured from the acquisition images (Section 2.2). 
No periodic signal is detected in the continuum variations. 
We note that, as the system became brighter, the central 
spike became harder to detect and the lines became slightly 
broader. The central spike is completely absent from the 
Hei A 6678 line on the third night (see Figure [3]), and barely 
detected in Ha. We interpret this behaviour as variability 
originating in the accretion disc, outshining the coronal or 
chromospheric emission from the white dwarf as it brightens, 
causing the central spike to disappear from the spectrum. 



6 For completeness, we note that the Ha emission line in the CV 
SDSS J003941. 06+005427.5 is triple-peaked ( |Southworth et al.| 
2010), but unlike the central spike observed in AMCVn stars, 
it has a large radial velocity amplitude (202 kms -1 ), which does 
not match the expected velocities of the white dwarf, the sec- 
ondary star, or the accretion disc bright-spot. jSouthworth et al.| 
(2010) suggest that it may originate from a coronal loop on a 
magnetically active secondary star, similar to those that have 
been observed in outburst spectra of IP Peg and SS Cyg jSteeghs 
et al.|1996| ). 
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3.3 Orbital period 

To measure the radial velocity variation, we used the dou- 
ble Gaussian technique of |Schneider &; Young| (p. 980) and 
Shafter (1983), as implemented in MOLL\|^] This measure- 



ment is sensitive to the high velocity line wings which orig- 
inate in the inner disc and is often used as an indication 



of the radial velocity of the white dwarf (e.g. Thorstensen 
[2000] ). 

We used two Gaussian functions of full width half max- 
imum FWHM = 400 kms -1 to account for the broad emis- 
sion lines, and varied the separation between the peaks from 
800 to 3000 kms -1 . For each resulting radial velocity curve 



we calculated the Lomb-Scargle periodogram flLomb|[l97^ 



Scargle 1982) and then folded the radial velocities on the 
strongest period. The folded radial velocity curve was then 
fitted with a circular orbit of the form 



V(t) = K sin [2tt(* - BJD )/P or b] + 7 , 



(1) 



where the reference time BJD is defined by the blue to red 
crossing of the measured velocities. 



A diagnostic diagram analysis (Shafter et al. 1986) 
showed that a/K is a minimum for a separation of 
1200 kms -1 , where a is the scatter around the fit. This sep- 
aration also gives a consistent orbital period and reference 
time from all four emission lines. We show the Lomb-Scargle 
periodogram calculated from the Hei A 7065 radial velocities 
in Figure [4] The strongest signal occurs at a frequency of 
22.075 cycles d -1 , corresponding to a period of 65.233 min. 
The nearest competing peaks are the ±1 day aliases at 
21.058 and 23.064 cycles d -1 . We carried out 10 000 boot- 
strap selections of the radial velocity curve, calculated the 
periodogram of each of these randomly selected subsets and 
recorded the strongest peak. 98.95 per cent of the radial ve- 
locity subsets resulted in 22.074 cycles d -1 as the strongest 
signal. 0.6 per cent of the subsets found the 21.058 cycles d -1 
alias marginally stronger, and only 0.45 per cent preferred 
the alias at 23.064 cycles d -1 . Thus we are confident that 
the strongest signal is the correct identification of the spec- 
troscopic period of this CV. The Hei A 7065 radial velocity 
curve, folded on this period, is shown in the middle panel of 
Figure [4] The parameters of the best fit of Equation [I] as 
derived from the three strongest lines, are listed in Table [3] 

The bottom panel of Figure [4] shows the Ha radial ve- 
locities folded on the same period. The amplitude of the 
variations is much smaller than those derived for the helium 
lines and the curve is highly non-sinusoidal. Its saw-tooth 
shape suggests that the emission region is non-circular. A 



7 molly was written by T.R. Marsh and is available from 
http : / / www. Warwick, ac . uk/go/ tr marsh / software / 
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Figure 3. Top row: Average, continuum normalised emission line profiles from the first (solid line) and third (dashed line) nights of 
observations on VLT/FORS2. The double-peaked profile is characteristic of Doppler motion in an accretion disc. On the first and second 
night, the Ha and Hel A 6678 lines additionally show a narrow component at the centre of the line. The lines also get broader as they 
get brighter, indicating that the emission originates from a higher velocity region. This is true even if the lines are normalised to the 
line peaks. The vertical axis is the same for all four line profiles, for ease of comparison. Second row: Trailed spectra from all three 
nights, normalised to the continuum and folded on the orbital period. Two orbital cycles are shown for clarity. The panels were scaled 
individually to the minimum and maximum flux in each panel. Third row: Doppler maps calculated from the first night's spectra. The 
central spike of the emission line maps to a bright spot at the origin of the map, consistent with the expected velocity position of the 
white dwarf. Bottom row: Doppler maps calculated from the third night's spectra. The white dwarf emission component has disappeared, 
as it is dominated by the disc emission. The Doppler maps are subject to an unknown rotation compared to the coordinate system such 
maps are usually displayed in (see text), as the zero phase of the white dwarf is not known. We used BJDq = 2455646.6655 as in Fig. [4] 



non-circular emission region is also supported by the shape 
of the Ha Doppler map (Section 3.5) in Figure [5] Even 
though we derive a consistent period from the four avail- 
able emission lines, the amplitudes of the variation differ by 
several tens of kms" 1 and in all cases are much larger than 
expected for the white dwarf in such a compact binary. This 
is a clear indication that the phase-dependent modulation 
of the line wings is not purely a result of the radial velocity 
of the white dwarf. The high velocity wings of the emission 
lines contain flux from other bright regions of the disc, e.g. 
the bright-spot, and cannot be used to track the motion of 
the white dwarf in this case. 

In an effort to isolate the bright-spot emission to mea- 
sure its velocity amplitude more accurately, we subtract 
from each spectrum the average spectrum from that night. 
This removes the constant flux, and leaves only the variable 
part of the emission. The signal-to-noise ratio is not high 



enough to fit these residuals in the individual spectra with 
a moving Gaussian to measure the period, but when folded 
on the period determined above, the disc bright-spot emis- 
sion clearly shows up in the Doppler maps and as an S-wave 
in the trailed spectra. The S-wave is also clearly visible in 
the continuum normalised spectra shown in the second row 
of Figure [3] Incidentally, folding the residuals on the either 
of the 1 day aliases, smears the emission out completely, 
so that no structure is visible in the trailed spectra. This 
is further confirmation that we have selected the correct 
alias as the orbital period. Measuring the bright-spot veloc- 
ity from the phase-binned spectra yields a semi-amplitude 
of 550 ±27 kms" 1 for the Hel A7065 line, 575 ±28 kms" 1 for 
the Hel A6678 line and 506 ± 16 kms" 1 for Ha. (Table^. 
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Table 3. Orbital properties as derived from each of the strongest emission lines. 



Ha Hei A 6678 Hei A 7065 



P orb (d) 0.04529(2) 0.04532(1) 0.04530(1) 

BJD 2455646.6660(3) 2455646.6663(2) 2455646.6655(2) 

K (km/s) 69.4 ± 2.9 285.8 ± 5.8 208.4 ± 5.1 

7 (km/s) 2.8 ±2.0 32.6 ± 4.7 16.1 ± 3.6 

Bright-spot velocity (km/s) 505.8 ± 15.7 575.0 ± 28.3 549.9 db 26.5 
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Figure 4. Top panel: Scargle periodogram calculated from the 
Hei A 7065 radial velocities of CSS1122-1110. The adopted orbital 
frequency is indicated by an arrow. Middle panel: The measured 
Hei A 7065 radial velocities folded on a period of 65.233 minutes. 
The dashed line is the best-fit sinusoid to the radial velocities 
and its parameters are shown in Table [3] Two phases are plotted 
for clarity Bottom panel: The Ha radial velocities, folded on the 
same period. Note the smaller amplitude compared to the He 
radial velocity variation. 



3.4 Binary mass ratio 

Superhumps are characteristic of dwarf novae in superout- 
burst ( |Whitehurst]|1988||Hirose fe Osaki||1990| ). A resonant 
interaction between the outer accretion disc and the donor 
star causes periodic flexing of the disc, which compresses 
the disc on the side opposite to the donor star and increases 
the viscous dissipation in this region. The observed bright- 
ness modulation is known as superhumps. During the later 
stages of the outburst, variable dissipation from the bright 
spot also contributes to the superhump signal, as it sweeps 
around the edge of the asymmetric disc. See Wood et al. 
(2011 ) for a detailed discussion of the origin of superhumps. 



The dynamical nature of this phenomenon means that there 
is good reason to expect the superhump period, P sh , to be 
related to the mass ratio q of the binary. Specifically, it is 



■-Porb)/-Porb, which is found 
Patterson et al.||2005|). This 



the superhump excess, e =(P S h 
to correlate well with q (e.g. 
empirical relation is derived from eclipsing dwarf novae, for 
which the mass ratio and superhump excess can be measured 
independently. 

There are three similar formulations of the e — q rela- 
tion which are commonly used, but unfortunately the low-e 
end of the relation is poorly constrained by observa tions. 
Patterson et al.| ( |2005| find e = 0.18g + 0.29<? 2 and |Kato| 
et al. ( |2009 ) give an updated version as e = 0.1 6g + 0.25g 2 . 
Both of these formulations assume that e = when q = 0. 
This is a reasonable assumption, since for a negligible sec- 
ondary mass, we would expect the tidal interaction with the 
disc to become negligible as well. Observations of the only 



known eclipsing AM C Vn system, SDSS J0 926+3624 flCop- 
perwheat et al.||2011), agree best with the Patterson et al. 



(2005) relation, so we favour that version in the calculations 
below. For completeness, we also note the third formula- 
tion of the e — q relationship, as given by Knigge] ([2006), 
q = (0.114 ± 0.005) + (3.97 ± 0.41) x (e - 0.025). Here the 
assumption e = = q is not made, so the value of q is 
probably overestimated for this binary. 

Using our spectroscopically determined period and the 



superhump period measured by Kato et al. (2010), we find 
a small superhump excess in CSS1122-1110, e = 0.0031 =b 
0.0007. This implies an extreme mass ratio, q = 0.017 =b 
0.004. For the average white dwarf mass found in CVs, Mi = 
(0.83 ± O.23)M ( |Zorotovic et al.|2011| ), this implies a very 
low mass donor, M 2 = (0.014 ± 0.005)M o . 



3.5 Doppler tomography 

The second row of Figure[3]shows the trailed spectra, centred 
on each of the emission lines covered by our FORS2 spectra. 
The spectra from all three nights are included in these trails, 
each normalised to the continuum and then folded on the 
orbital period. The central spike is visible as a faint line at 
the line centres of Ha and Hei A 6678, and an S-wave is seen 
in all four trails. 

We computed corresponding Doppler tomograms 
flMarsh & Horne||1988|) for each of the emission lines, for 



each night separately, as well as phase-binned together over 
the three nights. The increase in the line width as the system 
gets brigher blurs the map when the spectra from all three 
nights are included in a single map, so instead we show the 
maps calculated from the first and the third nights' spectra 
separately, in the third and fourth row of panels in Figure [3] 
The maps from the second night's spectra are near- identical 
to those from the first night. All maps were calculated using 
BJDo = 2455646.6655 as a phase zero reference time, as de- 
rived from the radial velocities of the Hei A 7065 line wings 



in Section [3.3| The resulting maps are clearly rotated with 
respect to the standard orientation usually used in these 
displays, which has the white dwarf and donor star aligned 
along the vertical Vx = axis ( |Marsh &; Horne||1988| ). In 
that orientation the disc bright-spot appears in the upper 
left quadrant of the map. This rotation seen here is a result 
of the unknown phase shift between our choice of B JDo and 
true zero phase of the white dwarf. The maps all show a 
similar structure, with a bright ring corresponding to emis- 
sion from the accretion disc, a bright-spot on the disc rim 
(equivalent to the S-wave in the trailed spectra) , and a small 
spot at the origin of the map, which corresponds to the cen- 
tral spike seen in the average Ha and Hei A 6678 line pro- 
files. We measured the velocity position of the central spike 
by fitting it with a two-dimensional Gaussian function. Its 
velocity is consistent with zero, so we are not able to use 
it to determine the zero phase orientation of the map. We 
have observed more than one binary orbit per night, so we 
have also calculated Doppler maps for each orbit individ- 
ually. From the position of the central spot in these, and 
the combined maps, we can constrain the radial velocity of 
the white dwarf to K\ < 16 kms -1 . On the third night of 
our observations, when the system was at its brightest, the 
emission was dominated by the accretion disc. As a result, 
the central spot disappeared completely from the Hei A 6678 
map, and only some weak emission remained at this position 
in the Ha map. 

As an illustration, we overplot the Hei A 6678 map with 
the Roche lobe parameters of a g = 0.017 binary in Figure [5] 
The assumed velocity positions of the white dwarf, donor 
star and velocity streams are also shown. We stress that, 
since the zero phase of the white dwarf is not known, we 
have the freedom to rotate the map about its origin, so this 
combination of parameters is not a unique description of this 
Doppler map. It only serves as an illustration of the extreme 
mass ratio q = 0.017, as determined from the e — q relation 
in the previous section. 



4 DISCUSSION 

Three channels are thought to contribute to the formation 
of AM CVn stars, (a) A double-degenerate channel, where a 
white dwarf accretes from another, helium-rich white dwarf 
(Paczyhski 1967; Faulkner et al. 1972). This ultra-compact 
configuration requires a second common envelope phase dur- 
ing the evolution to shorten the period to the observed val- 
ues, in addition to angular momentum loss by gravitational 
radiation. Potential direct progenitors for this type of binary 
are the detached double white dwarf systems, in particu- 
lar those with extremely low mass white dwarf companions 
( Brown et al.||20lTj 2012). (b) An alternative outcome of a 
double common envelope evolution is that the donor star 
is not fully degenerate when it leaves the common envelope 



(Savonije et al. 1986 Iben fe Tutukov 1987). The system 



starts mass transfer with the donor in the helium burning 
stage and the donor becomes increasingly degenerate as it 
evolves towards the period minimum of ~ 10 minutes, (c) 
If the initial separation is too large for a second common 
envelope to form, the system may evolve as a CV. An ul- 
tracompact binary may form if the initial mass ratio is such 
that the donor star undergoes significant nuclear evolution 
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Figure 5. Hei A 6678 Doppler map overplotted with the Roche 
lobes of a q = 0.017 binary. The assumed velocities of the white 
dwarf and donor star are labeled with x symbols and the centre 
of mass velocity at (0, 0) with a + symbol. The map was rotated 
by applying a —0.25 phase shift compared to the maps shown in 
Fig|3] 



before the onset of mass transfer ( Podsiadlowski et al. 2003) 



The density of the evolved donor allows the system to reach 
a much shorter orbital period than is possible for the hydro- 
gen CVs, with the exact value depending on the level of hy- 
drogen depletion in the donor when the mass transfer starts 



(e.g. Sienkiewicz 1984). The donor stars in these systems will 



have some hydrogen left in their envelopes as they evolve 
towards their period minima, but will become increasingly 
helium-rich during the evolution. This is the only formation 
channel in which the donor is predicted to have residual hy- 
drogen in its envelope, so the detection of hydrogen in the 
spectrum of an AM CVn binary conclusively identifies it as 
having evolved along this track ( |Pod siadlow ski et al.| 2003 
|Nelemans et al.||2010 ). 



The evolved CV formation channel is often dismissed as 
improbable or unimportant compared to the other channels, 
since it would require an evolutionary timescale comparable 
to, or even exceeding, the Hubble time to reach periods as 
short as those observed in AM CVn systems, if angular mo- 
mentum loss is driven by magnetic braking and gravitational 
radiation ( |Nelemans et al.||2010 Yungelson et al., in prep). 
From population synthesis models, |Nelemans et al.| |2004 ) 
estimate that less than two per cent of AM CVn binaries can 
form this way. However, if the initial mass ratio is q > 1, 
the binary will undergo a period of thermal timescale mass 
transfer (TTMT) until the mass ratio decreases to q < 5/6, 



allowing stable mass transfer (iTutukov et al. II 19851 Pylyser 



fe Savonije|1989| |Schenker et al.|2002| ). The rapid mass loss 
strips the donor from its hydrogen envelope, exposing the 
partially evolved core. When the system resumes normal 
mass transfer and angular momentum loss, it will appear 
essentially as a normal CV, except that the donor will be 
evolved compared to other CVs at the same orbital period. 
Those systems in which most of the hydrogen is lost during 
the period of TTMT, will be able to evolve to periods well 
below Pmin- Since the rapid mass transfer also removes angu- 
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lar momentum at a high rate, it is possible for these systems 
to evolve to ultrashort periods within a Hubble time. The 
fraction of CVs with evolved secondaries is expected to be 



high. Podsiadlowski et al. ( 2003| ) put this figure at ten per 
cent, and Schenker et al. ( 2002| suggest a fraction as high 
as a third. 

Observations of V485 Cen (59 min, Augusteijn et al. 
" 1996} and EIPsc (64 min, |Thorstensen et al.||2002[ 



1993 



support the idea that binaries which have undergone TTMT 
can evolve to an ultracompact configuration. The optical 
spectrum of EIPsc shows a donor star with spectral type 
K4±2 ( Thorstensen et al.| |20Q2), which is anomalously hot 
for a main sequence star at such a short orbital period. Its 
UV spectrum also reveals a low carbon to nitrogen abun- 
dance (Gansicke et al. 2003), indicating that the accreted 
material is CNO processed, so the white dwarf is accret- 
ing from the exposed, evolved core. Apart from V485Cen, 
EIPsc and CSS 1122- 11 10, there are two more short pe- 
riod binaries which are candidates for AM CVn progeni- 
tors of the 'evolved CV channel: CSS090331:102843-081927 



(52 min, |Woudt et aLl|2012| herea fter CSS1028-0819) and 
CSS111019:233313-155744 (62 min, |Woudt fc Warner|20lT 
hereafter CSS2333-1557). At the moment, little is known 
about either of these systems. CSS1028-0819 was discovered 
in outburst by the CRTS and subsequent outburst photom- 
etry revealed a period excess of e = 0.0530 flWoudt et al. 
20121 . VSNET observers ( |Kato fe Kinugasa|2009) reported 
the outburst spectrum to show both hydrogen and helium 
lines. No spectroscopic observations of CSS2333-1557 are 
available yet to determine whether it is an AM CVn system, 
a hydrogen accreting binary or a halo CV0 

There are only 42 semi-detached binaries known with 
periods below P^iL'- 36 AM CVn systems, 1 low metallicity 
halo CV, 3 ultracompact hydrogen accreting binaries and 
2 binaries discovered as a result of outburst activity, but 
still of unknown nature. The three confirmed ultracompact 
hydrogen accretors (CSS1122-1110, V485 Cen and EIPsc) 
therefore account for seven per cent of this total. If CSS 1028- 
0819 and CSS2333-1557 are confirmed as hydrogen accret- 
ing systems as well, this figure will rise to 12 per cent. Al- 
though these figures do not take selection effects into ac- 
count, such a high fraction of possible progenitors suggests 
that the evolved CV channel is non-negligible when consid- 
ering AM CVn formation. A statistically complete sample of 
ultracompact accreting binaries is not yet available. 

We show the cumulative period distribution of all ac- 
creting white dwarf binaries with known orbital periods in 
Figure [6] In the period range 59 — 66 minutes both long 
period AM CVn stars and short period hydrogen accreting 
systems are found, but currently there are no known sys- 
tems with periods between 67 and 76 minutes. This is sur- 
prising, because if the critical parameter which determines 
the minimum period a system can reach is simply the level 



8 One more short period hydrogen-accreting binary is known, 
SDSS J150722.30+523039.8 (67 min, |Littlefair et aLpOu?] ). It is 
excluded from the discussion, as it has been shown to be a CV 
in the Galactic halo, with a low mass Population II donor star 
(Uthas et al. 2011). Its short period stems from the lower atmo- 
spheric opacity of such a low metallicity star, which means that 
it has a smaller radius and fills its Roche lobe at shorter orbital 
periods. 
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Figure 6. Cumulative period distribution of all AM CVn stars, 
CVs with P rb< 90 min and short period hydrogen accreting bi- 
naries with known orbital periods. AM CVn stars span the period 
range 5 — 65 minutes, CVs the range 76 minutes and longer and 
the ultracompact hydrogen accretors 59 — 65 minutes. There are 
no known accreting white dwarf binaries with an orbital period 



between 67 and 76 minutes. The periods were taken from |Solheim 
( 2010 L |Levitan et al. 1(201 lj,|Ritter fe Kolb|(|2003t,|Gansicke et al 
( 2009), | Augusteijn et al | (1996), |Thorstensen et al.| ( |2002| >, |Lit- 
|tlefair et al.| ( |2007| ) and jWoudt et al.| ( |2012| >. Individual periods 
are marked by ticks. The unclassified binaries CSS1028-0819 and 
CSS2333-1557 are included as CVs, but highlighted by a small 
horizontal bar. We use an abbreviated notation to indicate the 
periods of V485 Cen, EIPsc and CSS1122-1110. 



of hydrogen depletion when the mass transfer starts, there 
should be no reason not to find binaries in this period range 
(see e.g. Sienki ewicz||1984 ). 

Unlike for EI Psc and V485 Cen, the secondary star is 
not detected in our spectra of CSS1122-1110 (see Figure [TJ, 
which means that it must be of later spectral type than the 
K star seen in EIPsc. The argument for a very low mass, 



late type donor (Section 3.4) is further supported by the 
mass ratios of those binaries, which are less extreme than 
that of CSS 1122-1110. Fo r V485Cen, we use the P sh mea- 
surement of Olech (1997) to calculate q = 0.13 =b 0.01, and 



for EIPsc, the P sh estimate from Uemura et al. (2002) gives 
q = 0.14 ± 0.01. The mass ratio we find for CSS1122-1110 
is much more extreme, q = 0.017 ± 0.004. It is however sim- 
ilar to the mass ratios of the long period AM CVn stars, 
e.g. GPCom (P orb =46.6 min, q = 0.018) and V396 Hya 
(P or b=65.1 min, q = 0.0125) (both taken from Nelemans 
2005), which suggests that it may be similarly evolved. Us- 



ing a mass-radius relation for a fully degenerate helium star 



(Eggleton 1986, |Verbunt fe Rappaport| [l988), the mass of 
a such a donor star at an orbital period of 65.233 minutes 
will be O.OO6335M , which translates to q min = 0.0076. This 
is even lower than the q measured for CSS 11 22- 11 10, in- 
dicating that the donor star in CSS1 122-1 110 is probably 
semi-degenerate . 

The observed properties of CSS 1122- 11 10 strongly sug- 
gest that it is an AM CVn system forming via the evolved 
CV channel. All but the most extreme of AM CVn stars 
which evolve along the evolved CV channel are expected 
to have hydrogen left in the donor star. The system is ex- 
pected to become increasingly helium-rich as it evolves to 
shorter periods, but a small amount of hydrogen could still 
be left in the donor envelope as it passes the period min- 
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imum. The strong (compared to AM CVn stars) hydrogen 
lines in the spectrum of CSS 11 22- 11 10 may suggest that is 
still evolving towards its minimum orbital period. The fact 
that it shows outbursts also supports this view, since the 
mass transfer rate must be high enough to trigger the disc 
thermal instability for outbursts to occur. This is in contrast 
with GP Com and V396 Hya which have never been observed 
in outburst, and are believed to have very low mass trans- 
fer rates. AMCVns showing outbursts are generally found 
in the 20-40 minute period range (e.g. Solheim 2010). On 



the other hand, the very low mass of the donor star could 
be interpreted as CSS 11 22- 11 10 having already evolved past 



its period minimum ( jPodsiadlowski et al.|2003[ Table Al). 

A detailed abundance study may help to determine 
the evolutionary status of this binary. The relative abun- 
dances depend on a number of factors, such as the hydrogen 
abundance at the onset of mass transfer, the accretion rate 
and the temperature structure of the disc. Element abun- 



dances are difficult to measure, but Nagel et al. (2009) suc- 



cessfully reproduced the spectral features of V396 Hya us- 
ing a non local thermodynamic equilibrium accretion disc 
model. The line ratios we measure from the SDSS spec- 
trum of CSS1 122-1 110 already suggests that the hydrogen 
content is very low compared to canonical CVs. We mea- 
sure a flux ratio Hei A5876/Ha= 1.49 =b 0.04, compared to 
Hel A5876/Ha:~ 0.2 — 0.4 for a typical dwarf novae (e.g. 
|Szkody|198ll|Thorstensen fc Taylor 1 200 1) . A more detailed 
abundance analysis requires a high signal-to-noise spectrum 
covering a wide spectral range, but a suitable spectrum is 
not yet available for this system. The relative strengths of 
CNO species in such a spectrum will further help to dis- 
tinguish between the different evolutionary scenarios (see 
Nelemans et al.|2010 ). 



5 CONCLUSION 

Time-resolved optical spectroscopy of the dwarf nova 
CSS100603:112253-111037 (CSS1122-1110) reveal an orbital 
period of 65.233 =b 0.015 minutes, well below the CV period 
minimum. As is usual for CVs, the spectra are dominated 
by Hei and Balmer emission lines, but the Hei lines are un- 
usually strong, likely a sign of high helium abundance. The 
Ha and Hei A 6678 lines display a narrow emission com- 
ponent, stationary at the line centre to within 16 kins" 1 . 
Using Doppler maps, we can associate this 'central spike' 
with emission from the white dwarf. Such a triple-peaked 
emission line structure, with a central spike originating from 
near the white dwarf, is a characteristic only seen in helium- 
dominated systems, not in CVs. We combine our spectro- 
scopically measured orbital period with the superhump pe- 



riod measured by Kato et al. (2010) to find an extreme mass 
ratio for this binary, q = 0.017 ± 0.004. Using the average 
CV white dwarf mass, this implies a very low mass donor 
star, M 2 = (0.014 ±0.005)M©. Comparing the properties of 
this ultracompact accreting binary with CVs, AM CVn sys- 
tems and hydrogen binaries below the CV period minimum, 
we argue that CSS 1122- 11 10 is the first compelling example 
of an AM CVn system forming via the evolved CV channel. 
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